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I. INTRODUCTION 

The Lovelock gravity [])] is a natural generalization of general relativity (Einstein gravity 
theory) in higher dimensional spacetimes. Its action is a sum of some dimensionally extended 
Euler densities, where there are no more than second order derivatives with respect to metric 
in equations of motion. Also it is known that the Lovelock gravity is free of ghost. Over the 
past years, due to the development of some theories in higher dimensional (> 4) spacetimes 
such as string theories, brane world scenarios etc, the Lovelock gravity has attracted a lot 
of attention, and some new features and properties have been revealed. 

The action of Lovelock gravity can be written as 

p 

£ = Y,CiCi, (1.1) 

where p < [{n — 1)/2]([A^] denotes the integer part of the number A^), q are arbitrary 
constants with dimension of [Length] n is the spacetime dimension and Ci are the Euler 
densities 

where the generalized delta function is totally antisymmetric in both sets of indices. From 
the Lagrangian (11. ip . one can get the equations of motion, Qab = 0, where 

P 1 

na raai - aifei---bi pcirfi jDCidi /i o\ 

^6 ~ 2«+l aifei aibi- y^-'^) 

1=0 

After introducing matter into the theory, we can get the equations of motion including the 
energy-momentum tensor of matter 

gab = 87TGTat. (1.4) 

From (11.31) . it can be shown that the equations of motion do not have more than second 
order derivatives with respect to metric and that the Lovelock theory has the same degrees 
of freedom as the ordinary Einstein gravity theory. Just because of this fact, the Lovelock 
gravity is free of ghost when expanded on a flat spacetime, avoiding any problem with 
unitarity [21]. Note that, Cq denotes the unity, and Ci gives us the usual curvature scalar 
term, while C2 is just the Gauss-Bonnet term. Usually in order for the Einstein gravity to 
be recovered in the low energy limit, the constant Cq should be identified as the cosmological 
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constant up to a constant and Ci should be positive (for simplicity one may take C\ = 1). 
For example, in n = 4, we have 



^CoQab + Ci{Rab - ^Rgab) 



8nGT 



ah 



:i.5) 



This is just the Einstein equation with the cosmological constant A = — Co/2 if we set Ci to 
be one. When n = 5 and p = 2, after expanding the Kroneker-Delta, one gets the equations 
of motion for Gauss-Bonnet gravity in the usual manner. 

In the literature on the Lovelock gravity, the most extensively studied theory is the so- 
called Einstein-Gauss-Bonnet (EGB) gravity. The EGB gravity is a special case of Lovelock's 
theory of gravitation, whose Lagrangian just contains the first three terms in (11. II) . The 
Gauss-Bonnet term naturally appears in the low energy effective action of heterotic sting 
theory ^ . Spherically symmetric black hole solutions in the Gauss-Bonnet gravity have been 



bund and discussed in 



IB, 3, 



and topological nontrivial black holes have been studied in 



3]. Rotating Gauss-Bonnet black holes have been discussed in js]. Some other extensions 
such as including the perturbative AdS black hole solutions in gravity theories with second 
order curvature corrections could be seen in ^]. In addition, the references in 10|] have 
investigated the holographic properties associated with the Gauss-Bonnet theory. And the 



papers m 



12 



13 



14l | gave some exact solutions for Vaidya-like solution in the Einstein- 



Gauss-Bonnet gravity. 

For a generic case, although the Lagrangian (II. ip looks complicated, some exact black 



hole so 



16 



13, 



utions have been found and their associated thermodynamics was investigated in 15|, 



18 



Q, 2^ . In the so-called third-order Lovelock gravity, that is, containing the first 
four terms in (11. ID . some exact solutions have been found in 21|. Furthermore, it is also 
known that those higher derivative terms in the Lagrangian (II. ip with positive coefficients 
arise as higher order corrections in superstring theories, and their cosmological implication 
has been studied 22 1 . 



In this paper, we are mainly interested in dynamical black hole solutions in Lovelock 



gravity by generalizing those discussions in 



11 



I2,ll3|,ll|. The organization of the paper 



is as follows. In Sec. II, we present a kind of generalized Vaidya spacetime in the general 
Lovelock gravity (II. ip . In Sec. Ill, we study thermodynamics of apparent horizon of the 
generalized Vaidya spacetime, and give corresponding entropy expression associated with 
the apparent horizon. In addition, we propose a way to obtain the generalized Misner-Sharp 
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energy in the Lovelock gravity. Section. IV. is devoted to conclusion and discussion. 



II. GENERALIZED VAIDYA SOLUTION IN LOVELOCK THEORY 

In the four-dimensional general relativity, the Vaidya spacetime is a typical dynamical 
one. In this spacetime there exists a pure radiation matter. The metric of the spacetime 
can be written as 

ds"^ = -f{r, v)dv'^ + 2dvdr + r'^d^l . (2.1) 

where / = 1 — 2m(f)/r, and dVL\ is the line element of a two-dimensional unit sphere. In this 
spacetime, the apparent horizon is given by / = 0, or r = 2m{v). The energy-momentum 
tensor for the radiation matter in the spacetime is given by Tab = fJ'^ah, where /„ = (!, 0, 0, 0) 
in coordinates (f , r, 9, 0). The quantity fi is the energy density of the radiation matter. Now, 
in an n-dimensional spacetime, assume a similar form of spherically symmetric metric 

ds'^ = -f{r, v)dv^ + 2dvdr + r^dnl_^ . (2.2) 

The energy-momentum tensor of radiation matter has a similar form as that in four dimen- 
sions. For this metric, it is not hard to calculate the nonvanishing components of Riemann 
tensor given by 

£/> £1 £1 

TDvr £ TDvi __J sri r)ri __J sri 

2 ' 2r ' ~ 2r ^ ' 

ri ft vi i- ^ - f i- 

^"vj = -^^j ' ^"ri = , R'\i = —^^ki ■ (2-3) 

Here a prime/overdot denotes the derivative with respect to r /v. Substituting these results 
into the equations of motion (11. 4p . and using identities 

and 

^a,---a^.ra^^b ,b ^ 2[n - (m - 1)1 [n - (m - 2)]5"^■■■;'"-^ 



we find the equations of motion 
(n-2)! 



Ci 



'n-2i-l)\ 



0, (2.4) 
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J_ 

2r 



l-f 



i-l 



(2.5) 



and 



i 

(n-3)! 



1] (n-3)! fl-f 



+ Ai- 



2 I (n- 22 -3)! V 



n-2i-2)\\ 2r \ r 



+ 2i- 



(n-3)! 



n - 2z - 1)! V 2 J \ r 



f"\ fl-f 



+ Ai{i - 1) 



(n-3)! 



{n-2i- l)\ \ 2r 



f \ fl-f 



i-2 



(2.6) 



Other components of are given by = Qy, Q^! = 0. The components Qj are not 
independent, because they are hnearly expressed in terms of drQ^ and Q^, as will be shown 
below. Defining a new function F{v,r) 

1 - f{v,r) 



F(v, r) 



(2.7) 



we can put the equations (12. 4p and (12. 6p into the forms 

^ 2 n-2i-l !r"-2 L J 

and 

- V c ^ Tr"-! " - 



(2.8) 



(2.9) 



^From these two equations, it is easy to show ^] = 5] [rdrG^^ /{n — 2) + Q"!^] , so ^] = do not 
yield independent equations. Integrating the equation (12.81) leads to an order-p algebraic 
equation for F{v,r) or f{v,r) 



E 



{n-2)\ ■ _ lQ'KGm{v) 



(2.10) 



(ra - 2i - 1)!^ l^„-2r"-i ' 

where m{v) is an arbitrary function of v, which appears as an integration constant. (Cer- 
tainly, to ensure some energy condition, this mass function should be positive.) The coeffi- 
cient lQ7[G/Qn-2 is chosen such that m can be interpreted as the mass of the solution when 
m is a constant. Using 

tF'-^F = t\-^\{^^] , (2.11) 



d 



dv 
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in the equation (12. Sp . we have 



8nGf^ = yc. }''''^^-\, ^F\ (2.12) 
^ ^ 2{n-2i-l)\dv ^ ' 



Comparing equations (12.101) and (12.121) . we obtain 



miv] 



Thus, we obtain a kind of radiating Vaidya spacetime in a generic Lovelock gravity by solving 
equation f l2.10l) . When the dimension of spacetime is five, that is, n = 5, the solution reduces 

nnnn 

to the one reported by some authors in references [HI, 112, 113|, HJ] , while in = 4, the solution 
is just the familiar Vaidya solution of general relativity. 

Now we further generalize the Vaidya spacetime in Lovelock gravity to more general case. 
Note that for the metric (12.21) we have — S^, so the energy-momentum tensor of matter 
has to satisfy = T^. Certainly, the matter of pure radiation discussed above satisfies the 
constraint. In fact, they are = = 0. If we further assume the spherical part of the 
energy- momentum tensor has the form T- = crT^ = aT^ (where o" is a constant, and the 
repeat index i does not sum), then from the equation VaTj^ = or the explicit expressions 
of ^5 in equations (12.41) . (12.51) and (12.61) . we can find 

n — 2 

d,T: + drTi + r; = , (2.14) 



and 



a.T; + ^!i^^Ki^T; = o. (2.15) 

r 



As a result, for the pure radiation matter with = = 0, one can find that has to be 
proportional to l/r"~^. This is consistent with the equation (I2.13p . Next, in the case with 
= 0, the equation (I2.15P tells us that and have the form 

T; = T; = C(t;)r^("-2)(^-"\ (2.16) 

where C{v) is a function of v. The off-diagonal part of the energy-momentum tensor T^", 
i.e., the component has to satisfy the equation ( I2.14p . Now the equations of motion 
= = 87rGC(t;)r-("-2)(i-a) j^odify the equation to 

^^^- 2(n-2/-l)! t^""^^]' = C(.)r("-)^ (2.17) 
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Integrating this equation, we have 

^ (n-2i-l)! V^n-2r"-i ^n-i y ' v ; 

where m(f ) is an arbitrary function, appearing as an integration constant again. Here, 
= J drA^~'^^", which is 

e(r)=ln(r), (2.19) 

for cr = — l/(n — 2) and 

(n-2)(T+l 

otherwise. Certainly, to ensure some energy condition for the energy-momentum tensor, the 
parameter a and function m(t;) and C{v) should satisfy certain consistency relations. For 
example, C{v) < and — 1 < cr < 0. These relations have been discussed in One can 
see from the relation (12. lip that TJ' is given by the partial derivative of the equation (I2.18p 
with respect to v. Thus we have 



_ ~ _ mjv) , C{v)e{r) 



This is consistent with the equation (12.141) . So the energy-momentum tensor of matter in 
this case can be written as 

Tab = P^Lh - PiLn-b + nJb) + aPqab , (2.22) 

where Ua is a null vector which satisfies lati"" = — 1. In coordinates {f,r, ■■■}, we have 
/a = (1, 0, 0, ■ ■ ■ ), and = (//2, —1, 0, ■ ■ ■ ). The tensor qab is a projection operator which is 
given by qab = Qab + la^b + ^b^a- So the metric (12. 2p can be put into the form Qab = hab + qab-, 
where 

hab = -Lnb - hua (2.23) 

is the metric of the two-dimensional spacetime transverse to the (n — 2) -dimensional sphere. 
Certainly, in coordinates {f , r, ■ ■ ■ }, the line element of hab can be expressed as —f{v, r)dv'^ + 
2dvdr. The quantity P is the radial pressure with the form P = —C{v)r^^'"^'^^^^^'^\ 

Generally, it is not easy to solve equation ( 12.18P when p is greater than one. However, 
for some special cases, one can explicitly get analytic solutions. Now, we give some simple 
examples. 
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1^1). Gauss-Bonnet theory: The simple radiating Vaidya solution in the Gauss-Bonnet 



gravity without a cosmological constant has been found in the references 



Here we can directly solve the second order algebraic equation fl2.18p . and give the solution 



'f{v, r)dv^ + 2dvdr + r^dnl^ 



2 5 



where f{v,r) satisfies the equation fl2.18p . which in this case becomes 



F + a{n-3){n- 4:)F^ 



IGnG 



n 



m[v) 



2 \^n-2r 



C{v)Q{r) 



(2.24) 



(2.25) 



Here, we have set Cq = 0, Ci = 1 and C2 = a. This parameter a, called the Gauss-Bonnet 
coefficient, has the dimension of length squared. Since it is an algebraic equation of order 
two, the solutions have two branches in general. Only one branch has the general relativity 
limit and is given by 

^2 



f{v,r) 



1 + 



2(n-3)(n-4)a 



1 - 




647r(n-3)(n-4)a / m{v) C{v)e{r) 



n 



(2.26) 

where, for simplicity, we have set G = 1. From f{v,r) = 0, we can obtain the trapping 
horizon, which is also its apparent horizon in this case. The radius of apparent horizon 
satisfies 



„n— 3 



r'^ " + a{n — 3){n — 4) 



n— 5 



16n 



n - 2 \a 



miv) 



civ)eirA: 



'n-2 



(2.27) 



Thus, in general, the radius of the apparent horizon is a function of v. This solution is the 
same as the one in [13!] if a cosmological constant is included. 

(2). Dimensionally continued Lovelock gravity: An interesting case is to choose some 

nnnn 

special values for the coefficients Cj, i = 0, ■ ■ ■ ,p. In [15|, 116|, 117|, 118| a set of special coefficients 
has been chosen so that the equation (12.181) has a simple expression. In odd dimensions the 
action is the Chern-Simons form for the AdS group, while in even dimensions it is called 
Born-Infeld theories constructed with the Lorentz part of the AdS curvature tensor. In the 
odd-dimensional case, say n = 2p-|- 1, i.e., Chern-Simons theory, we can choose 



(n-2i-l)!r ]£ 



-n+2i 



(2.28) 



where the parameter i is a length scale. Note that Cj differ from the coefficients in the 
references [15, 



16 



17 



as 



181] only by a factor of {n — 2i)\. At the same time we choose l/167rG 



(n-2)! 



(2.29) 
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Then the equation (12.181) gives the solution 

2 

fiv, r) = 1 - [m{v) + ^]„„2C(^;)e(r)]^; + ^ , 
For the even- dimensional case, say = 2p + 2, we set Cj, z = 1, ■ ■ ■ , p, to be 

Ci = {n — 2i — 1) 
and the gravity coupling constant 



-n+2i 



11-2 



levrG (n-2)!' 
Then the equation (12 .181) gives the solution 



f{v,r) 



1 - 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



The solutions (I2.30p and (I2.33P reduce to the static cases if m{v) and Civ) do not depend on 



coordinate v. In these cases, the solutions have a unique AdS vacuum [l5|, [ly, [l7|, [l8| . These 
solutions fl2.30p a nd (12.330 with vanishing C{v) have already been found by M. Nozawa and 



H. Maeda in 



(3). Pure Lovelock gravity: In this theory 



23 



24 



251], only two of coefficients Cj are 



non-vanishing: one is Cq and the other is Ck with 1 < k < p. We can normalize them as 
co/(n — l){n — 2) = — 1/£^ and Ckin — l)\/{n — 2k — 1)! = where I and a are two 

length scales. Then the equation (12.180 becomes 



a 



2k-2 Tpk 



pl^ = — 



levrG 



P n-2 \nn-2r 



miv) 



C{v)Q{r) 



r,ri~l 



(2.34) 



Thus, in the case of even dimensions, if the right hand of (12.341) is non-negative, we have 

.2 r^,2 ifi^^.,2 / ^^^^ C{v)Q{r 



f{v,r) 



1 ± 



+ — 



n 



^n-2r 



n—l 



+ 



r,n— 1 



while, in the other case of odd dimensions, 

„2 r^,2 i«_r^^2 



f{v,r) 



1 + 



r 



a" 



OL" IGnGa 

+ 

n 



m{v) C(f)0(r) 



^n-2r 



n-l 



+ 



(2.35) 



(2.36) 



Thus, for = 1, one gets a generalized Vaidya black hole with a cosmo logical constant. 
When m{v) and C{v) are two constants, these solutions reduce to the static solutions. 

For dynamical black holes, it is difficult to study their thermodynamical properties. In 
four- dimensional Einstein gravity, Hayward has proposed a method to discuss this issue 30, 



31 



32^, 



331]. In the next section, we will discuss the thermodynamics of the apparent horizon 



of the new solutions in this section. 
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III. ENTROPY AND ENERGY OF THE APPARENT HORIZON 

In the early 1970s, it was found that four laws of black hole mechanics in general relativity 
are very analogous to four laws of thermodynamics. Due to Hawking's discovery that black 
hole radiates thermal radiation, it turns out that it is not just an analog, but they are 
indeed identical with each other. Then thermodynamics of black hole has been established 
soundly (for a review, see j26|). In black hole thermodynamics, the temperature and entropy 
of a black hole are given by Teh = K/27r and S'eh = ^/4, where k and A are the surface 
gravity and area of event horizon of the black hole, respectively. In higher derivative theory 
of gravity, such as Gauss-Bonnet gravity, it turns out that the area formula for black hole 
entropy no longer holds. In fact, black hole entropy gets corrections from these higher 
derivative terms. For a diffeomorphism invariant theory, Wald [^^j showed that the entropy 
of black hole is a kind of Noether charge, and obtained a formula for black hole entropy, 
now called Wald formula, associated with the event horizon of black hole. For more details. 



see the review paper by Wald 26 1 



In the literature, most of discussions on thermodynamics of black holes have been focused 
on stationary black holes. Thermodynamics of a black hole is associated with the event 
horizon of the black hole, which is the boundary of the past of future infinity. Therefore the 
event horizon is a null hypersurface and depends on some global structures of the spacetime, 
so it is difficult to study the event horizon of a dynamical (time-dependent) spacetime. In 
general relativity, there exists another kind of horizon, named apparent horizon, which is not 
heavily dependent of global structure of spacetime. In the standard definition of apparent 
horizon, one first slices the spacetime, and then finds the boundary of trapped region in 



each slice. This boundary (two-dimensional surface) is called the apparent horizon 28|. 
Also the three-dimensional hypersurface of the union of all these two-dimensional surfaces is 
usually called the apparent horizon. Over the past years, several generalizations of horizon 
have been proposed, such as the trapping horizon by Hayward, the isolated horizon and the 
dynamical horizon by Ashtekar et al. The relations and differences among those horizons 
have been discussed in a recent review {23]. For a dynamical black hole, the outer trapping 
horizon and the dynamical horizon are not null hypersurfaces but spacelike hypersurfaces of 



the spacetime 



291]. Therefore, some well-known results associated with the event horizon. 



such as Wald entropy formula, may not be applicable to those horizons. 
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In this section, we will discuss the entropy and the energy associated with the apparent 



horizon of the generalized Vaidya solution in Lovelock gravity 



theory; 



(I2.18p . Before going 



on, let us first give a brief review on the work of Hayward [30|, iMl, 13^, l33|] . Although his work 
focuses on four-dimensional Einste in g ravity, it can be straightforwardly generalized to higher 
dimensional Einstein gravity 3J, ISSj. These discussions reveal the deep relation between 



equations o 



also 



36 



37|, 



motion and thermodynamics of the spacetimes. For relevant discussions, see 



38 



39, 



40 



4l|. For an ra-dimensional spherically symmetric spacetime {Ai,gab) 



we can write the metric in the double null form 



ds'^ = habdx°'dx^ + r^(x)(ifi^„2 i 



(3.1) 



where {x"} are coordinates of the two-dimensional spacetime (M, hab) which is transverse to 
the (n — 2)-dimensional sphere, r(x) is the radius of the sphere. Similarly, one can also divide 
the energy-momentum tensor of matter into two parts. One part denoted by Tab (do not 
be confused with the total energy-momentum tensor) corresponds to the two-dimensional 
spacetime hob- From this energy- momentum tensor, one can define two important physical 
quantities: the work density W = —l/2h°'^Tab, which corresponds to the work term in the 
first law, and the energy supply \I'a = T^dbr + Wdar. By using these two quantities and the 
Misner-Sharp energy inside the sphere with radius r, 

(n-2), 



E 



IQtv 



-f^n-2r"-^ (1 - h-'dard^r) 



(3.2) 



one can put some components of the Einstein equations into the so-called unified first law 



dE = + WdV . 



(3.3) 



where A = Vtn-2^"'~'^ and V = Qn-2f"'~^/{n — 1) are the area and the volume of an {n — 2)- 
sphere with radius r. For a vector ^ tangent to the trapping horizon of the spacetime, 
Hayward showed that on trapping horizon, one has 32 1 



A^I/aC = —C.A = —C.S = TSS . 



(3.4) 



where S = A/A, T = h/2tt, and k, is the surface gravity defined by k = DaD"'r/2. Here, Da 
is the covariant derivative associated with metric hab- The 6 operator in the right hand side 
of the equation (13.41) should be understood as follows: take a Lie derivative with respect 
to ^, and then evaluate it on the apparent horizon. Note that the left hand side of the 
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equation (13 ■4p represents the amount of energy crossing the trapping horizon. Therefore the 
equation (13. 4p may be understood as the Clausius relation for mechanics of dynamic black 
holes: 5Q = T6S with 6Q = A'^a^"'- By projecting the unified first law onto the trapping 
horizon, the first law of thermodynamics is given by 

6E = T6S + W6V . (3.5) 

In the generalized Vaidya spacetime found in the previous section, the apparent horizon 
is just a kind of trapping horizon, so we do not emphasize the difference between these two 
concepts in the present paper. For the generalized Vaidya spacetime (12. 2p . the apparent 
horizon is given by f{v,r) = 0. From the definition of surface gravity, it is easy to find 
K, = /'(f,ryi)/2. We now discuss thermodynamics on apparent horizon/trapping horizon 
for the dynamical solutions in Lovelock gravity theory. Following [^], we may move all 
terms of Qab except for the Einstein tensor into the right hand side of the field equations 
Qab = —^T^Tab and rewrite them in the standard form for Einstein gravity 

Gab = SnTab = Svr (t^™^ + T^?) . (3.6) 

Here, T^^™'' is the energy- momentum tensor of matter (12.22^ . The effective energy-momentum 
tensor T^l' has the form 

8nT^l^ = Hab = gab-Gab, (3.7) 

where Gab is the Einstein tensor of the spacetime. In this way, we can go along the line 
of Einstein gravity, although we are discussing a gravity theory beyond the Einstein grav- 
ity. Thus, the work term and the energy supply defined on the two-dimensional spacetime 
{M^,hab) are, respectively, 

W = -^h'^'^Tab = IV^"") + Vr^'^ , 

■^a = T%r + Wdar = + , (3.8) 

where W^"^^ and \E'i"'-* are defined by using T^f, while ly(^) and ¥a^ by T^^^'\ and the 
contraction is taken over the two-dimensional spacetime hab- It is easy to find in our case, 

W^"^^ = -P = C(t;)r"("-2)(i-<x) ^ ^3_g^ 
so the energy supply for the matter is given by 

v^i-)=/2/,. (3.10) 
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Now, because we are treating an effective Einstein gravity (i.e. view the non-Einstein term 
Hab as an energy-momentum tensor of higher curvature terms), the Misner-Sharp energy 
(13. 2p is apphcable. For the generahzed Vaidya spacetime, we can replace the term h°'^dardbr 
in equation (13.21) by /(f , r). It is easy to see that equations (13. 8p and (13.21) satisfy the unified 
first law (13.31) . 

On the apparent horizon, we have the Clausius-like equation 

A^ai" = + A¥^^C = ■ (3.11) 

Stt 

where ^ is a vector field tangent to the apparent horizon. This vector can be determined 
as follows. Since / = on the {n — 2)-dimensional marginal surface, the Lie derivative of 
/ with respect to ^ should vanish, so £5/ = on this surface. This means that on the 
apparent horizon we have 

Cdj + Cdrf = 0. (3.12) 

Noting that n = f'{v,rA)/2, we have = —f /2k,. In the case of / = d^f = 0, we have 

= 0. Therefore on the apparent horizon, it is reasonable to set C,a to be 

Here, the normalization of C, is not important and will not play a key role in the following 
discussion. It is easy to find C,°'C,a = — //^- So the generator or the tangent vector of the 
apparent horizon is not a null vector unless / = 0, i.e., the static case. On the other hand, 
for the static cases, one has = ra" = (d/dv)"" on the apparent horizon, as expected. 

Note that the heat flow 6Q is determined by the pure matter energy-momentum tensor or 



35|. Thus, on the apparent 



pure matter energy-supply. This point has been emphasized in 
horizon we can rewrite the equation (13.111) as 

6Q = A^"^^™) = ^6 A - AC¥f . (3.13) 

An interesting question is whether the right hand side of the equation (13.131 ) can be cast into 
a form T6S of the Clausius relation as in the Einstein gravity. The answer is affirmative, 
since the right hand side of the (I3.13P indeed can be written in a form of T6S, where 
5* will be given below. One can get this entropy expression by directly substituting the 



expression of Hab in equation (13.70 into the equation (13.130 . In fact, the authors in [35 1 
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obtain an entropy expression of apparent horizon in Lovelock gravity, but in a setting of 
FRW universe. Here let us stress that it is not always possible to rewrite the right hand side 
of the equation (13.131) in a form T6S, scalar-tensor theory and f{R) gravitational theory 
being the counterexamples 35 1. 

In order to find the entropy expression on the apparent horizon, here we use a little 



35( 1 . By using the explicit form of in equation (13.101) , 



different method from the one in 
we have 

6Q = = AjllaC = -Aj2 . (3.14) 

/,From equations (EHH), fl2:T8D and (ESID, it is not hard to find 

/ Svr^ (n-2i-l)! ^ ' 

Hence we have on the apparent horizon, 

SQ = -Af^ = ^(l](^y c. /^"''^, r^-^A = ^C,S = ^6S , (3.16) 
^ ^ 27r I 2k M 4 ^ (n - 2z - 1)! / 27r^ 27r ' ^ ^ 

where all the calculations are done on the apparent horizon. This is very similar to the 
Clausius relation 6Q = T6S if we define the temperature by T = K/27r and the entropy of 
apparent horizon 5* by 

Ay p-2)l 

The entropy (13.171) is the same as that in the static cases if one replaces ta with the event 
horizon radius of static black holes in Lovelock gravity 19|. 

Now, we have shown that A'^^"^^ term can be written as T5S term in the first law. Can 
the equations of motion in Lovelock gravity be rewritten as the unified first law (13. 3p in 
Einstein gravity theory, where the left-hand side is completely determined by spacetime 
geometry, while the right hand side is determined by matter in spacetime? We can rewrite 
fl33D as 

dE - A"^^^^ - W^^UV = v4^(™) + W'^^'^dV . (3.18) 

Note that the left hand side of the equation is totally determined by geometry because 
^('^^ and ly(^) are defined by geometric quantities as well as E. This implies that the left 
hand side of the equation (13.181) can give us an energy form like the Misner- Sharp energy in 
Einstein gravity. Indeed it is not hard to show that the left hand side of the equation can 
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be cast into a form of dEi, where the function El is given by 

i ' 

Thus in the Vaidya-hke spacetime for Lovelock gravity theory we arrive at a generahzed 
unified first law 

dE^ = + W^'^UV . (3.20) 

After projecting on the apparent horizon, we obtain the first law of apparent horizon 

6El = T6S - P6V . (3.21) 

where we have used W^"^^ = —P. If we write the function f(v, r) back into the form 
h°'''dardbr, the energy fl3.19p is nothing but the generalized Minser-Sharp energy suggested 



43 



44j . where the generalized Misner- Sharp 



by H. Maeda and M. Nozawa in recent papers 
energy is shown to be a quasilocal conserved charge associated with a locally conserved 
current constructed from the generalized Kodama vector. Here let us stress that our proce- 
dure to find the entropy and the energy expressions is a direct generalization of the method 
proposed by Hayward to Lovelock gravity theory. The two methods, though different, lead 
to the same result. It would be interesting to use our method to study the unified first law 
in other gravity theories and to get the corresponding Miner-Sharp energy. 

In some sense. Lovelock gravity theory is special since this theory is ghost-free, and the 
equations of motion are in fact second order for metric, as in the case of Einstein gravity. 
These properties may help us separating matter from geometry effectively. For other theories 
without such properties, such as f{R) gravity, the relation between equations of motion and 
thermodynamics is complicated. Recently, Eling, Guedens and Jacobson have pointed out 
that equations of motion for f{R) gravity correspond to a non-equilibrium thermodynamics 
of spacetime. To get correct equations of motion, an entropy production term has to be 



added to the Clausiu relation |37l]. However, to what extent can the spacetime be described 
by equilibrium thermodynamics is still an open question. 

In summary, the definition of heat variation 6Q is important in our procedure. The 
above discussion shows that 6Q should be defined by pure matter energy supply or pure 
matter energy-momentum tensor, through which we can find a reasonable entropy associated 
with the apparent horizon. With this entropy and the unified first law in Einstein gravity 
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theory, we can obtain a proper energy form, i.e. generalized Misner-Sharp energy, completely 
determined by spacetime geometry. Thus we can establish the generalized unified first law 
for Lovelock gravity theory. 

The generalized unified first law (13.211) of the apparent horizon is a version of physical 
process for the first law of black hole thermodynamics, i.e., an active version of the first law. 
In this version of the first law, one considers a spacetime and matter energy-momentum 
tensor, and calculates the variation of thermodynamical quantities which are induced by the 
matter (it enters into horizon). The first law then establishes a relation among the variations 
of those thermodynamic quantities. In the passive version or phase space version of first law, 
on the other hand, one compares the thermodynamical quantities of two spacetimes differing 
from each other by a small amount of variables, and then gives variations of these quantities, 
which lead to the first law. This physical process of the first law can shed some light on 
the passive version of the first law. For the solution (12.181) . here we give some suggestion 
on the passive version of the first law. By using equations ( 12.18^ and the definitions of the 
apparent horizon and the surface gravity, the surface gravity can be written as a function 
of the radius of the apparent horizon 

1 



^(n -2i- l)c,r^2^+2 ^ i6vrC(t;)ri"-')'^-"+' 



(3.22) 



where ta = ta{v) and the coefficients q are defined by 



^.r^^^- (3.23) 



(n-2i- 1)! 

The energy inside the apparent horizon is given by 

M, = ^j2c^rl~''-\ (3.24) 

i 

and the pressure on the apparent horizon by 

P = -C(t;)r-("-'^('"") . (3.25) 

Then it is easy to check that the following relation holds 

dMi = TdS - PdV . (3.26) 

where the entropy is given by (13.171) . while dr^i denotes a variation of the apparent horizon, 

i.e., from ta to ta + dr^- This is a variation of the state parameter ta of the dynamical 



17 



black hole. It should be noted that the variation operator "d" is not an exterior derivative 
"d" of the spacetime manifold, but a difference of two nearby points in the solution space 
of the theory. Equation fl3.26p is just the passive version of the first law. Another point we 
wish to stress is that the physical meaning of energy Ml is ambiguous. It is not an ADM 
mass of the spacetime (although it is an ADM mass in the static limit), but a generalized 
Misner-Sharp energy inside the apparent horizon. By using equation (I2.18p . the energy 
can be expressed as 

Ml = m{v) + ^]„_2C(^;)e(r^) . (3.27) 

Thus, we see that the generalized Misner-sharp energy inside the apparent horizon for the 
solution f l2.18p gets the physical meaning as a generalized Bondi mass if the matter is pure 
radiation. However, for other type of matters, such as the term leading to the nonvanishing 
C(f), the Bondi mass should be modified. Further, it turns out difficult to define the Bondi 
mass at null infinity in odd dimensions. This has been pointed out in references Q, I47I . 
For the static case, ta = t+, equations (13.220 and (13.260 reduce to the first law of black hole 



thermodynamics in Lovelock gravity theory 



19|. 



IV. CONCLUSION AND DISCUSSION 

In this paper, we have obtained a kind of generalized Vaidya solutions (I2.18P in a generic 
Lovelock gravity. Explicit forms of these solutions are given for some special cases, such as 
Gauss-Bonnet gravity, dimensional continued gravity and pure Lovelock gravity. We have 
also investigated the unified first law on the apparent horizon of this spacetime, and given 
the entropy (I3.17P associated with the apparent horizon and the energy (I3.19P for (n — 2)- 
dimensional sphere with radius r. In our procedure, we treated all terms of Lovelock gravity 
except for the Einstein tensor in equations of motion as an effective energy-momentum tensor 
in Einstein gravity theory, to which we could apply the unified first law proposed by Hayward 
for Einstein gravity theory. Defining the heat fiux by pure matter energy-momentum tensor, 
we obtained an entropy expression (I3.17P associated with the apparent horizon. Substituting 
this entropy into the unified first law leads to an energy form (I3.19P for Lovelock gravity. 



which is the same as the one suggested by H. Maeda and M. Nozawa [43|, |4J]. This method 
can be further extended to other theories of gravity, for example, brane world theory. 

Here a comment not discussed in section HI is in order. We did not discuss whether the 
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apparent horizon (or trapping horizon) is "outer" and "inner" [s^. In general, for higher 
dimensional cases, the parameter p in (11.11) is more than one, so the equation (I2.18P is a 
higher order algebraic equation of f{v, r) and may have multi-horizons. We have assumed 
in section III that the apparent horizon is an outer apparent horizon (trapping horizon) of 
black hole spacetime with positive surface gravity l/2/'(f , r)|r=T-^ > 0. In fact, the result 
of this paper also holds for the cosmological event horizon (inner apparent horizon) with 
negative surface gravity l/2f'{v,r)\r=r^ < 0. In this case the Hawking temperature of the 
apparent horizon is T = |fi;/27r|. 

Note that, while investigating the thermodynamical properties of the new solution, we just 
focussed on the mechanics of the dynamical black hole. To establish the thermodynamics of 
apparent horizon, one has to show that there exists Hawking radiation with the temperature 
T = k/27c as in the case of static black holes. In fact, this may be proved by 



recent works on Hawkin g ra diation and entropy of dynamical spacetime 
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bl 
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owmg some 
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51 



52| 



(for a recent review, see |53l|). 

In a dynamical spacetime the apparent horizon is not necessarily the same as the event 
horizon. For the solution (I2.18p . the apparent horizon is simply given by an algebraic 
equation /(f,r^(f)) = 0, from which the event horizon is very different. The event horizon 
requires the global information of the spacetime. However, if we assume the event horizon 
is a null hypersurface of the spacetime, it is not hard to find that the event horizon of 
solution ( I2.18P is given by /(f,r+(f)) = 2dr^{v)/dv. Therefore, in general, it is not easy 
to give an explicit expression for the location of event horizon because we have to solve 
a differential equation instead of an algebraic equation. For some simple cases, such as 
the four-dimensional Vaidya spacetime and the five- dimensional Vaidya spacetime with the 
Gauss-Bonnet term, we may give the explicit formula of the location of the event horizon 
for some simple mass functions m{v). Taking the four-dimensional Vaidya spacetime in the 
Einstein gravity as an example, the event horizon is given by r_|_ = 2m(f)/(l — 2r+). Some 
authors have calculated the associated temperature and entropy by using brick wall model 



and gravity anomalies method 5J, l55| : 



S = — = nrl . (4.1) 



87im(v) 47rr+ 4 

At this stage a serious problem may arise. For a dynamical spacetime such as the Vaidya 
solution, which horizon, the apparent horizon or the event horizon or both of them, is indeed 
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associated with Hawking radiation and the entropy? This issue has been commented in 56 1 
(references therein). But it is fair to say that this is still an open question. Obviously, this 
is a quite interesting issue worth further studying. 
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